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Abstract. Recently spatial as well as temporal variations of the fine structure 
constant a have been reported. We show that a "runaway domain wall", which 
arises for the scalar field potential without minima, can account for such variations 
simultaneously. The time variation is induced by a runaway potential and the spatial 
variation is induced by the formation of a domain wall. The model is consistent with 
the current cosmological data and can be tested by the future experiments to test the 
equivalence principle. 
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1. Introduction 

String theory is the most promising approach to unify all the fundamental forces in 
nature. It is believed that in string theory all the coupling constants and parameters 
(except the string tension) in nature are derived quantities and are determined by the 
vacuum expectation values of the dilaton and moduli. However, only few mechanisms 
(see for example, pQ) are known how and when to fix the dilaton/moduli. On the 
other hand, we know that the Universe is expanding. Then it is no wonder to imagine 
the possibility of the variation of the constants of nature during the evolution of the 
Universe. 

In fact, it is argued that the effective potentials of dilaton or moduli induced by 
nonperturbative effects may exhibit runaway structure; they asymptote zero for the 
weak coupling limit where dilaton becomes minus infinity or internal radius becomes 
infinity and symmetries are restored in the limit [21 [3] . Thus it is expected that as these 
fields vary, the natural "constants" may change in time and moreover the violation of 
the weak equivalence principle may be induced [3], SI El El E] - 

Hence, any detection or nondetection of such variations at various cosmological 
epochs could provide useful information about the physics beyond the standard model. 
In this respect, the recent claims of the detection of the time variation [U |9l [10] as 
well as the spatial variation [TT] of the fine structure constant a may hint towards new 
physics. 

Narrow lines in quasar spectra are produced by absorption of radiation in 
intervening clouds of gas, many of which are enriched with heavy elements. Because 
quasar spectra contain doublet absorption lines at a number of redshifts, it is possible 
to check for time variation in a simply by looking for changes in the doublet separation 
of alkaline-type ions with one outer electron as a function of redshift. 

Webb et al. [8] introduced a new technique (called many-multiplet (MM) method) 
that compares the absorption wavelengths of magnesium and iron atoms in the same 
absorbing cloud, which is far more sensitive than the alkaline-doublet method. From the 
latest analysis of Keck/HIRES (High Resolution Echelle Spectrometer) 143 absorption 
systems for 0.2 < z < 3.7, they found that a was smaller in the past [TO] : 
Act 

= (-0.543 ±0.116) x 10" 5 . (1) 

a 

Moreover, recently, Webb et al. analyzed a dataset from the ESO Very Large Telescope 
(VLT) and found the opposite trend: a was larger in the past [UJ. Combined with the 
Keck samples, they claimed the spatial variation of a [TTJ : 

= (1.10 ± 0.25) x 10~ 6 (r/Glyr) cos 9, (2) 

a 

where r is the look-back time r = ct(z) and 9 is the angle between the direction of the 
measurement and the axis of best-fit dipole. 

However, concerning the claimed time variation of a, similar observations from 
VLT/UVES (Ultraviolet and Visual Echelle Spectrograph) have not been able to 
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duplicate these results [T21 EH dU US HE]. It is to be noted, however, that the analysis 
by Srianand et al. [12] may suffer from several flaws. For example, the uncertainty 
in wavelength calibration in [12] may not be consistent with the error in Aa/a JTT] . 
According to the analysis of the fundamental noise limitation [17], the systematic errors 
in [12] may be several times underestimated. Recent detailed re-analysis of Srianand et 
al. and Chand et al. confirms these concern: flawed parameter estimation methods in 
a x 2 minimization analysis [18] (see however, [19]) and systematic errors in the UVES 
wavelength calibration [20] . 

Overall, although the claims of the detection of the variations are not confirmed by 
independent methods (however see [21]), the claims are also not disputed seriously. We 
regard the claims of the spatial/temporal variations currently refuse to deny or confirm. 
In this paper, we take the claims of the spatial/temporal variations of a seriously and 
attempt to explain the data qualitatively. 

One major difficulty in explaining both temporal (Eq. ([1])) and spatial (Eq. fl2])) 
variations is that the spatial variation across the horizon scale (Aa/a ~ 1CT 5 at 
r ~ 14Glyr) is as large as the time variation during the Hubble time. If the time variation 
is induced by a scalar field [221 E31 EH [25j [26], such a cosmologically time evolving scalar 
field is very light with its mass being comparable to the Hubble parameter and hence its 
relative fluctuation is very small. So it is almost impossible to explain simultaneously 
both temporal and spatial variations by the light scalar field and its fluctuation. One 
possible solution to this problem is to consider nonlinear objects like topological defects 
[27] or giant voids [28]. In this paper, we consider the former possibility since the scalar 
field is dispersive and does not trace the matter density perturbation muchj|| Ref. [30] 
considered a domain wall to explain the spatial variations only but did not consider the 
time variations. We point out a certain type of domain walls can be utilized to explain 
not only the time variations of a but also the spatial variations of a of the same order 
of magnitude as the time variations. 

The paper is organized as follows: In Sec. 2, we detail our model and then in Sec. 
3 we study several constraints on the model parameters. Sec. 4 is devoted to summary. 



2. Runaway Domain Wall and Varying a 



For definiteness, we consider the theory described by the following action 



S 



d 4 xy/—g 



pi 



R 



- 2 mf 



V(<f) 



167rct 



£?(0) V 



Here M p i = 1/ VSttG = 2.4 x 10 18 GeV is the reduced Planck mass, <j) is the real symmetry 
breaking field, ao is the bare fine structure constant and S m denotes the action of 
other matter (relativistic/non-relativistic particles and dark energy). We consider the 

| It is possible to induce the spatial variations by the environmental dependence [24j [29], but it is 
difficult to explain the variation of dipole type like Eq. ([2]). 
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following scalar field potential of runaway type: 

M 2p+i 

Even if the potential has no minimum, the discrete symmetry <H- — can be broken 
dynamically, which results in the formation of a domain wall. Such defects are dubbed 
"vacuumless" defects in [3T], but we prefer to call them "runaway" defects since they 
arise from the potential of runaway type. 

For 3> a, from the balance between the kinetic energy and the potential energy, 
one finds — (M p+2 x) 1 ^ p+1 ^ and the energy density is proportional to x~ 2p ^ +1 \ where 
x is the distance from the wall [31]. We assume p > 1 so that the tension of the wall 
is finite. Then, the width of the wall is estimated as 5 ~ o~ p+1 /M p+2 and the tension of 
the wall is given by 

M 2P+4 M p+2 

a ~ - — 6 ~ -. (5) 

r a 2 P a p-i v i 

The profile of such a runaway domain wall solution is shown in Fig. [TJ In the 
cosmological situation, we replace x with the Hubble distance if -1 so that 

oc a 3 W p+ ^ (6) 

during the matter dominated era in accord with the tracker solution for V oc 0~ 2p [32] . 
This scaling solution is useful to account for the cosmological time variation of a; in 
the opposite case of <fi < a, the scalar field near the local maximum exhibits thawing 
behavior [33] and moves very slowly, and hence it is difficult to explain the cosmological 
time variation of a. 

The coupling function B(<p) in front of the electromagnetic kinetic term induces 
the spatio-temporal variations of a since a(<f>) = ao/B(<p). We consider the following 
coupling function : 

5(0) = e-^pi. (7) 
Since the effective a for small £ is given by 



oM^l + f^j, (8) 
the time variation in either side of the wall is given by 

- = ^ = ± , 3 (9) 

where we have used Eq. (jfj]), and the spatial variation across the wall is given by 
Act A0 

Thus, the opposite time variation of a between the Keck (increasing a) and the VLT 
(decreasing a) as well as the spatial variation of the same order of magnitude as the 
time variation during the Hubble time (Aa/«o ~ \a\faQH~ 1 ) are accommodated in this 
model. 
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Figure 1. <j) (upper) and the energy density (lower) of a runaway domain wall for 
p = 2 in a flat space. 



For definiteness, we shall choose p = 2 henceforth so that <\> oc a 1 ! 2 . Then the 
measured a by the QSO absorption lines ([8j |9l [101 [HJ [151 El HI]) are fitted as a 
function of a 1 / 2 as (a(a) — ao)/ao ~ —6.2 x 10 _6 (1 — a 1 / 2 ) and hence 



— -3.1 x KTV/ 2 #. 
Using this fit in Eq. (jHJ) at z ~ 2, we find that 



3 x 10" b . 



pi 



Putting this into Eq. fflO|) . it implies that 
Aa 



~ 7 x 10" 



(12) 



(13) 



which explains naturally the largeness of the spatial variation. 

Before closing this section, we would like to comment on the formation and 
dynamics of domain walls. For \<p\ a, the (tree-level) potential can be expanded 
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as 

V(cj>) ~ V - \m\<? + + ■ • • , (14) 

where V = M 8 /a A , ml = 4M 8 /a 6 , and A = 12M 8 /cr 8 . By taking into account the finite 
temperature effects, the effective potential around the origin reads 

V T (<f>) = V(<f>) + e^rf + \t^ 2 + ■ • • , (15) 

Thus, the phase transition occurs around the critical temperature T c ~ ^cr so that 
domain walls are formed. After some relaxation periodjjj domain walls evolve according 
to the linear scaling solution [31] , in which typical scale of the domain wall is comparable 
to the Hubble scale and its energy density is roughly given by 

p waU ~ fiH oc -. (16) 
3. Experimental Constraints 

Let us now discuss the observational and experimental constraints on the parameters. 



Sachs- Wolfe Effect: A domain wall induces the temperature anisotropy by the Sachs 
Wolfe effect [27]. The gravitational potential due to the wall at the horizon scale is 
27rG/!-£f ~ 1 ~ (l/4)M p J 2 /iiJ _1 which induces the temperature anisotropy via the Sachs- 
Wolfe effect. The requirement that this should be less than 10~ 5 gives 

M<30GeV (is) 1/4 - < 17 » 

The present energy density of a runaway domain wall within the horizon scale is esti- 
mated as [iHq 2 / (4:7tHq 3 /3), which should be much less than the critical density SM^Hq. 
This also gives a similar bound as Eq. (117jl . 



The Violation of the Weak Equivalence Principle: Since the effective mass is very 
light: W 77 ~ M 4 /0 3 < Kr 41 GeV(M/10GeV)V/10 15 GeV)- 3 for (j> > a, the scalar 
(f) mediates a long-range force via the coupling to nucleons, leading to the violation of 
the weak equivalence principle j22j [23]. The modification of the nucleon mass follows 
from the electromagnetic corrections. To leading order in a these corrections are given 
by [35], 

5m p = B p 5a/ao = 0.63MeV5a/a , (18) 

Sm n = B n 5a/ao = — 0.13MeV<5a/a , (19) 

§ The friction force due to the thermal plasma is estimated as rri^T 2 v (v : wall velocity) and is always 
subdominant in comparison to a force per unit area ~ a/t coming from the curvature ^ t. Thus, the 
friction effects on the domain wall dynamics due to the thermal plasma are negligible. This is simply 
because the field <f> consisting domain walls only weakly interacts with the thermal plasma. 
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where m p and m n are the proton and the neutron masses, and B p ~ 0.63 MeV and 
B n ~ —0.13 GeV are the Born terms for the proton and the neutron, respectively!]]] 
Hence, from Eq. (jSJ), the exchange of induces a composition dependent long-range 
force. A test body of mass m experiences the acceleration induced by the 0-exchange 
force [221 123] 

e 



AuM^mr 2 



{N^B P + N*B n ) {N p B p + N n B n ) , (20) 



where N pn (N P)7l ) are numbers of protons and neutrons in the Earth (the test 
body), in addition to the usual Newtonian acceleration due to the Earth mass Me' 
g = Me I '87T M 2 f 2 . The difference in accelerations between the two test bodies 



a 



in Eotvos-Dicke-Braginsky type experiments is parametrized by the Eotvos ratio: 
rj = 2\a\ — a.2 |/|«i + Q2 1) where a\ and a 2 are the accelerations of two bodies. Here 
we assume that the test bodies have almost equal masses m\ ~ m 2 , which implies 
N Pi i + N ns i ~ N Pj 2 + N n> 2. In the present case, rj is estimated as 

^_^,_ 2 e f N p E B p + N*B n \ / AN p B p + AN n B n \ 



a g m 2 y N r f + N* J \ N pA + N nA 

where we have used m, ~ {N n ^ + N p ^)m and Me — (N^ + N^m with m being the 
atomic mass unit (~ 0.931 GeV) and AN P = N Pi \ — N Pt2 , AN n = N n ^ — N nt2 . Adopting 
the typical values for N pn and AN PtTl [221 E3] , we find 



p,n 

2 



V* 3 x 10 " 14 (ib^ J • (22) 

This should be smaller than the current experimental bounds rj < 2 x 10~ 13 [36], which 
gives 

f < 2.6 x 10" 3 . (23) 



Allowed Parameter Region: The present value of is estimated as 

0o * (M 4 ^ 1 ) 1 / 3 . (24) 
From Eq. ([12]) and Eq. we obtain 

M ~ 30GeV 7 . (25) 

Moreover, the requirement of > cr to account for the time variation of a, from Eq. 
(IM]), leads to 

M>6GeV (Tra) 3/4 - < 26) 

For example, for a ~ 10 15 GeV, M ~ 30GeV and £ ~ 10" 3 satisfy Eqs. (TT71). (125) . 
and (126]) . This gives O — 8 x 10 15 GeV. In Fig. [2] we show the allowed range of the 
parameters M and a together with the relation Eq. ( 125]) . 



|| There also exists corrections to the binding energy of nucleon [35l [22] , but the effects do not change 
the result much. 
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Figure 2. Allowed parameter space. Upper region is excluded due to the Sachs- Wolfe 
effect (or large density parameter) Eq. (fT7|) ; lower region is excluded because of (j> < a 
and the absence of the scaling solution, Eq. (|26|) . Dotted lines explain the QSO data 
(Eq. ([25])) with ( = 2x 1(T 3 , 1(T 3 , 3 x 1(T 4 , 1CT 4 from bottom to top. 



Window j or Future Experiments: From Eq. f|T7|) and Eq. f[2"6"j) . upper bounds on M 
and a are found: 

M < 70GeV and a < 2 x 10 16 GeV, (27) 
which imply a lower bound on £ from Eq. ( I25jl : 

3 x 10~ 4 < £ < 2.6 x 10" 3 . (28) 
This in turn provides a window for r\ from Eq. ( l22l) : 

3 x 10^ 15 <?|<2x lir 13 . (29) 

Therefore, orders of magnitude improvements of the experimental limits on the weak 
equivalence principle by proposed experiments (such as MICROSCOPE [37], SR-POEM 
|38j . Galileo Galilei [39] and STEP [ID]) could lead to the detection of the violation 
of the weak equivalence principle induced by the scalar exchange force or refute this 
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model. These experiments (in particular, MICROSCOPE launched in 2012) can test 
the violation of the equivalence principle better than t] = 10 -15 . Therefore, the model 
can be tested within a few years by these gravitational experiments. 



4. Summary 

Motivated by possible detections of spatial and temporal variations of a, we have 
proposed a model based on a domain wall of runaway type. We have found that it 
is possible to construct a model to explain both variations simultaneously. We have 
studied the cosmological constraints on the model and found that the model can be 
made consistent with the current cosmological data and can be falsified by the future 
experiments to test the equivalence principle. We note that a model is not limited to a 
runaway potential, but we can construct a model with local minima so that the vacuum 
expectation value is determined by a runaway dilaton [4"T] . 

We have focused on a in this paper since our prime purpose was to provide a 
existence proof of a model. However, unless forbidden by symmetry, the direct couplings 
of 4> to fermionic matter should exist, which result in the violation of the equivalence 
principle and in spatio-temporal variations of the proton-to-electron mass ratio. There 
are some indications of a non-zero value of a spatial variation of it [32] • It would be 
interesting to study the consequences of such matter couplings further. 
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